Generic First-Order vs. Continuous Quantum Nucleation of Supersolidity 
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We analyze the nucleation of supersolid order out of the superfluid ground state of bosons on 
the triangular lattice. While the stability of supersolidity against phase separation in this system 
is by now well established for nearest-neighbor and long-range dipolar interactions, relevant for 
two-dimensional arrays of ultra-cold polar molecules, here we address directly the nature of the 
superfluid-to-supersolid transition. Based on symmetry arguments and quantum Monte Carlo sim- 
ulations, we conclude that this quantum phase transition is driven first-order beyond the line of 
particle-hole symmetry. Along this line, the transition is continuous and its scaling behavior consis- 
tent with the three-dimensional (3D) XY universality class. We relate this finding to a 3D 7L§ clock 
model description of the enlarged symmetry of the solid order parameter field. In the generic case 
however, the symmetry reduces to that of a 3D Z3 clock model, which reflects the first-order nature 
of the generic superfluid-to-supersolid quantum phase transition on the triangular lattice. 

PACS numbers: 67.80.kb, 75.40. Cx, 64.70. Tg, 75.40.Mg 



Polar molecules [l| are considered promising candi- 
dates for the realization of novel quantum states of mat- 
ter. They feature inherently long-ranged dipolar interac- 
tions, and state dressing by static electric and microwave 
fields allows to tune the inter-particle potential over a 
wide range and even into a regime that is dominated by 
three-body interactions [2, 3] . One feasible scenario is the 
realization of a supersolid phase upon loading bosonic 
polar molecules on a triangular optical lattice In 
fact, ensembles of hetero- nuclear molecules, such as KRb, 
RbCs, or LiCs in the rovibronic ground-state have been 
produced an d recently the superfluid to Mott in- 

sulator transition of 87 Rb atoms has been demonstrated 
on a triangular optical lattice [loj . Already for nearest- 
neighbor repulsive interactions, geometric frustration on 
the triangular lattice gives rise to an extended supersolid 
regime at sufficiently low kinetic energy 



3.5 




1st order 



T — ' — I — ' — T 



-1 — . I 1 L_ 



3.25 3.5 3.754 

Ou/V) c 

J 1 L 



0.1 0.12 0.14 0.16 0.18 

t/V 



0.2 



11-1 



This 

supersolid phase was found to be stable under extended 
dipolar interactions Q and for strong three-body repul- 
sions that emerge between dressed polar molecules 17 1. 

While such setups thus appear promising for an exper- 
imental realization of a supersolid state of matter, the 
nature of the transition from the superfluid to the su- 
persolid regime, at which solid order nucleates, has not 
been systematically addressed thus far in previous nu- 
merical studies. Initial works concluded the transition at 
or close to half-filling to be continuous, without however 
specifying the universality class of the quantum phase 
transition [H, [H, HH . O n the other hand, a recent clus- 
ter mean-field study concluded that the quantum phase 
transition is driven first-order beyond half-filling, and ex- 
hibits anomalous hysteresis behavior [l9j]. This calls for 
a careful re-examination of the quantum phase transition 
towards supersolidity on the triangular lattice. 

Here, we employ large-scale quantum Monte Carlo 
(QMC) simulations to assess the above mentioned sce- 
nario of a first-order superfluid-to-supersolid transition. 



FIG. 1. (Color online) Ground state phase diagram for hard- 
core bosons with nearest neighbor repulsion on the triangu- 
lar lattice near the superfluid-to-supersolid transition region. 
The inset shows the jump A of the order parameter corre- 
sponding to the solid order shown in the upper inset along 
the transition line as a function of the critical ratio (fj,/V) c . 
The dashed line in the inset is a quadratic fit for A < 0.05 to 
the deviation of (fJ,/V) c from 3 (half- filling). Errors are below 
the symbol size. 



We find that the transition is indeed first order beyond 
the particle-hole symmetric line of half-filling, while at 
half-filling the transition is continuous, and consistent 
with the 3D XY universality class. We relate our numer- 
ical findings to the peculiar properties of supersolidity on 
the triangular lattice, the corresponding discrete ground 
state degeneracies and the symmetry of the effective or- 
der parameter field, to describe the transition in terms of 
a g-state, Z g , clock model with q = 3 (6) for the generic 
(particle-hole symmetric) case, respectively. 

In the following, we consider the basic model of 
hard-core bosons on the triangular lattice with nearest 
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neighbor repulsions, described by the Bose-Hubbard like 
Hamiltonian 



(i) 



H = -t^2(b]bj + h.c.) + V ^2 riiUj - fj, 

(ij) (ij) i 

Here, hi (b\) annihilates (creates) a particle at lattice 
site i and n 2 ; = b\bi denotes the local density opera- 
tor. The tunneling matrix element t connects nearest- 
neighbor sites (ij), and V denotes the repulsion between 
adjacent sites. Furthermore, /i denotes the chemical po- 
tential controlling the filling n of the lattice. Before 
presenting our numerical results on the superfluid-to- 
supersolid transition, that are summarized in Fig. 1, we 
perform a symmetry analysis of the solid order across the 
superfluid-to-supersolid transition. The conclusions from 
these considerations are then shown to be consistent with 
the results from quantum Monte Carlo simulations. 

The phase diagram of the model in Eq. (1) has been 
subject to various analytical [ll|, EM 0] and numerical 
studies [l2 - 15 1 . For small values of t/V, two incompress- 
ible solids form, with \/3 x V3 diagonal long-range order 
(DLRO) of filling 1/3 and 2/3 respectively, related by 
particle-hole symmetry around fi/V — 3. For large val- 
ues of t/V, a uniform superfluid (SF) persists, where the 
system exhibits off-diagonal long-range order (ODLRO) 
and a finite superfluid density p s . Supersolid (SS) phases 
with both DLRO and ODLRO appear between the two 
solid lobes and prevail down to t/V = at p/V = 3. 
These supersolid phases are characterized by local den- 
sity fluctuations (61,62,63) = (=F2e, ±e', ±e') around a tri- 
angular plaquette, where gj = (rij — (n)) . The different 
signs relate to supersolids of filling n < 1/2 (SS-A) and 
71 > 1/2 (SS-B), respectively. At half filling (p/V = 3), 
both supersolid phases coexist at a first-order SS-SS tran- 
sition where both S/N and p s are continuous, whereas 
the pattern of the density fluctuations around a plaquette 
changes from SS-A to SS-B 0. 

In terms of the phase 9 of the complex order parameter 

ee 10 = ex + g 2 e~ i27r / 3 + g 3 e i2jr / 3 , 

the two supersolids are distinct in that (cos(30)) < 
for SS-A and (cos(30)) > for SS-B [lj]. Furthermore, 
this quantity respects the three-fold rotational symmetry 
of the lattice (9^9 + 2w/3). We thus expect that the 
nucleation of solid order at the generic (i.e. p/V 7^ 3) 
SF-SS quantum phase transition relates to that of a 3D 
three-state, Z3, clock model. The Z3 clock model equals 
a three-state Potts model 



21] 



which in 3D undergoes 
a first-order ordering transition |22j. Thus the generic 
SF-SS quantum phase transition is expected to be first- 
order. 

At p/V = 3 however, where the SS-A and SS-B co- 
exist, one finds (cos(30)) = 0, while (cos(60)) > 0. The 
symmetry is thus enhanced to a six-fold rotation of the 
order parameter by 7r/3. This symmetry enhancement 
has consequences on the nature of the SF-SS phase tran- 
sition, since the symmetry breaking now relates to that 
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FIG. 2. (Color online) Structure factor S/N for t/V = 0.15 
obtained at j3 = 100 for L — 12 and 36. Inset: histograms of 
S/N across the SF-SS transition for L = 36. 



of a six-state, Iiq, clock model. The 3D Z§ clock model 
has a continuous transition within the 3D XY universal- 
ity class [2TI. [23l425l| . Hence, we conclude that the SF-SS 
quantum phase transition at /i/V = 3 is continuous and 
within the 3D XY universality class. Thus far, in our 
analysis we neglected the coupling between the solid or- 
der and superfluidity; however, as argued in Ref. [ioT ]. 
the 3D XY fixed-point is stable in the presence of such 
interactions. From these results, we expect that the line 
of first-order SF-SS transition has a critical endpoint at 
fi/V = 3. In this scenario, the continuous transition 
at the particle-hole symmetric point of half-filling is not 
driven first-order via fluctuations [27j |. but turns discon- 
tinuous due to a reduced order parameter symmetry be- 
yond the line of half-filling. 

We next employ unbiased QMC simulations based on 
a dircctcd-loop algorithm in the stochastic series expan- 
sion representation (28l - [3Cl | to assess the SF-SS transi- 
tion. The interactions are decoupled into three-site tri- 
angular plaqucttcs, which allows for efficient updates. 
Our simulations are performed on systems with linear 
system size L up to 54 (the number of sites N = L 2 ), 
employing periodic boundary conditions in both lattice 
directions. Ground-state properties are obtained by tun- 
ing the temperature T sufficiently low, as detailed be- 
low. The superfluid density p s is obtained from mea- 
suring the winding number fluctuations 
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is detecting by measuring the corresponding structure 



factor S 



rj '{rii m) at wave vector 



H3 " v~> 3 ' 

Q = (27r/3,0), where the position of lattice site i is de- 
noted iv Since the problem is particle-hole symmetric, 
we need to consider only the regime p/V > 3. 

Distinguishing a (weakly) first-order from a continu- 
ous (quantum) phase transition can be challenging due 
to restrictions in the accessible lattice sizes with respect 
to the relevant correlation length. We thus start our 
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FIG. 3. (Color online) Filling n, superfluid density p s and 
structure factor S/N for t/V — 0.13 for various system sizes 
as functions of p/V across the SF-SS transition. 



FIG. 4. (Color online) Histograms of S/N obtained at t/V = 
0.13 for different system sizes. Inset: Extrapolation of the 
peak position A to the thermodynamic limit. 



analysis in a regime well away from half-filling, i.e. in 
the vicinity of the solid phase. Fig. 2 shows the structure 
factor S/N for t/V = 0.15, upon driving the chemical po- 
tential through the SF-SS quantum phase transition. A 
pronounced jump in (S/N) and p s (not shown) develops 
as one increases the system size, characteristic of a first- 
order transition. The inset in Fig. 2 shows histograms 
recorded in the close vicinity of the transition point, dis- 
playing a robust two-peak structure at (p/V) c ps 3.7545. 
This indicates that the ordered and disordered phase co- 
exist at this point. The jump in the order parameter of 
about A ps 0.04 is estimated from the position of the 
peak at finite S/N, cf. also the inset of Fig. 1. 

Moving in parameter space towards p/V = 3, the 
DLRO (and also ps) weakens, and the jump A in the or- 
der parameter decreases, as seen from the inset of Fig. 1. 
In particular, Fig. 3 shows the filling n, the superfluid 
density p s and the order parameter S/N as functions 
of fx/V for t/V = 0.13, in the vicinity of the SF-SS 
quantum phase transition. The data appears to vary 
smoothly across the transition, in pronounced contrast 
to the results shown in Fig. 2 for t/V = 0.15. The 
histograms of S/N Fig. 4 for sufficiently large systems, 
however, reveal the emergence of a two-peak structure 
at (p/V) c ~ 3.48075, which remains stable in the ther- 
modynamic limit, as shown in the inset of Fig. 4. We 
estimate the jump in the order parameter at the tran- 
sition point for t/V = 0.13 as A ps 0.0125 form Fig. 4. 
We also performed simulations at fixed t/V = 0.12 and 
fj,/V = 3.1, but we were not able to observe pronounced 
features in the S/N histograms apart from only a slight 
broadening of the peak in S/N, at least within the ac- 
cessible systems sizes up to L = 60. In fact, based on 
an approximate quadratic scaling A cx ((p/V) c — 3) 2 ob- 
served in the numerical data near half-filling (cf. the fit in 
the inset of Fig. 1), we conclude that in this regime the 



first order transition has become extremely weak, with 
an estimated A < 0.002 at t/V = 0.12, and is thus es- 
sentially resolution limited on the accessible system sizes. 
We thus obtained clear evidence for a first-order SF-SS 
transition over a wide range of the phase diagram away 
from half-filling, with a decreasing discontinuity upon ap- 
proaching the half-filled limit. Next, we focus on the 
SF-SS transition directly at p/V = 3, where based on 
symmetry considerations the quantum phase transition 
is expected to be continuous. Consistent with previous 
simulations 12-14j|, we obtain no indication for a first- 
order transition at p/V = 3 within the range of consid- 
ered system sizes (up to L = 54). In order to assess, if the 
transition indeed is of 3D XY type, we employ finite-size 
scaling analysis. The finite-size data of S/N for differ- 
ent system sizes in the vicinity of the SF-SS quantum 
phase transition is shown in the left inset of Fig. 5. Re- 
lated to the dynamical critical exponent z = 1 [26}, we 
performed the simulations fixing TL = 0.1V. We ex- 
tract the critical exponents using a conventional finite- 
size scaling ansatz for the structure factor S/N and the 
susceptibility X = (W) / dr(Ei,j e iQ ^- r ^ 71^)^(0)) 
at Q = (2tt/3,0), where n t (r) = e lHT fiie' lHT . The finite 
system data for S/N near a quantum critical point should 

follow the scaling relation S/N = L~ 2 ^"g (^~L l / v ^j . 
Similarly, we extract 7 and v from the scaling of %■ Here, 
t = t/V and t c denotes the position of the quantum crit- 
ical point. From our analysis, we locate the critical point 
at (t/V) c = 0.1108(2) and obtain the critical exponents 
v = 0.67(2), j3 = 0.32(2) from the scaling of the struc- 
ture factor, and (t/V) c = 0.1105(3) and v = 0.68(4), 
7 = —1.3(2) from the data collapse of the susceptibility. 
The error bars were obtained using standard bootstrap- 
ping, cf. the right inset of Fig. 5 for typical distributions. 
These results are consistent with recent high-precision 
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FIG. 5. (Color online) Data collapse of the structure fac- 
tor S/N (left panel) and the susceptibility \ (right panel) at 
fi/V = 3 for system sizes from L = 12 to 54 at LT = 0.1V . 
Left inset: finite-size data of S/N. Right inset: bootstrapping 
histogram for the critical exponents /3 and v from S/N data. 



values /3 = 0.3486(1), 7 = -1.3178(2) and v = 0.6717(1) 
of the 3D XY universality class [32j . The raw data and 
the data collapse using the reference values of the critical 
exponents are shown in Fig. 5. 

In conclusion, we analyzed the nucleation of superso- 
lidity of hard-core bosons on the triangular lattice; based 
on a symmetry analysis and quantum Monte Carlo simu- 
lations, wc found that this transition is first-order beyond 
the line of half-filling, whereas at half-filling, the transi- 
tion is continuous and consistent with the 3D XY uni- 
versality class. It will be important to revisit the quan- 
tum phase transitions between supcrfluid and supcrsolid 
phases also on other lattice geometries: For example, 
in case of a spin- 1/2 bilayer system, the superfluid-to- 
supersolid transition was reported from numerical simu- 
lations to belong to the 3D Ising universality class [33| , in 
contrast to the field-theoretical expectation for a checker- 
board supersolid state [3, 34 1. It will be particularly 



interesting to address the question, if the transition be- 
tween the checkerboard supersolid state and the super- 
fluid is in fact driven first-order by fluctuations for generic 
fillings, as suggested from an e-expansion [26l| . or if the 
exotic non-Bosc-liquid fixed-point obtained from a fixed- 
dimensional renormalization group analysis [26j could be 
realized. Our results imply that on the triangular lattice 
supersolid domains proliferate upon driving a supcrfluid 
system across the transition to the supcrsolid region - 
the need to distinguish such a phase coexistence from a 
mere mixture of (non-supcrfluid) solid and (non-solid) su- 
pcrfluid domains (like at a first-order supcrfluid-to-solid 
transition) could well challenge a robust experimental 
identification of the emerging supcrsolid state. 
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